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The bosonic excitations induced by the u) meson propagation in dense nuclear matter is studied 
within the framework of random phase approximation. The collective modes are then analyzed by 
O^l ■ finding the zeros of the relevant dielectric functions. Subsequently we present closed form analyt- 

ical expressions for the dispersion relations in different kinematical regime. Next, the analytical 
behaviour of the in-medium effective propagator for the uj meson is examined. This is exploited to 
calculate the full spectral function for the transverse (T) and longitudinal (L) mode of the uj meson. 
(3J[)' In addition, various sum rules are constructed for the uj meson spectral density in nuclear medium. 

Results are then discussed by calculating the residues at the poles and discontinuities across the 
cuts. 
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One of the cardinal challenges facing the nuclear physics today is to determine the properties hadrons at finite 
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temperature and/or chemical potential. Such studies are important both in the context of high energy heavy ion 
collisions (HIC) and fixed target electron scattering experiments |l],|2| . Medium effects have also been found to be 
important to explain pion induced pion production data. It is seen that a clear enhancement of tt + tt~ pair in the low 
invariant mass region could be accounted for by invoking the in- medium a meson spectral function Similarly 
it is observed that the photoabsorption data in 7 + A system can be explained by invoking the medium modified p 
O \ meson spectral function [EJ. 

In heavy ion collision, among others, the light vector mesons have acquired particular interest. This is because 
of their dileptonic decay channels. Dileptons, once produced leaves the medium without suffering further scattering 
^ . and hence provide a penetrating probe to study the properties of the dense nuclear medium as has been discussed 
J^j ■ at length in the recent literatures |pj. Experimentally the dilepton invariant mass spectra measured at CERN/SPS 
suggests that in-medium properties of the vector mesons, in particular that of p meson, are crucial to understand the 
■ low mass enhancement of the lepton pairs. 

In the theoretical arena such studies were stimulated by the pioneering work of Brown and Rho ^| where it was 
proposed that, in hot and/or dense matter, vector meson masses would drop from their values in free space as a 
precursor phenomenon of chiral symmetry restoration. Since then, several theoretical models have been invoked to 
study the in-medium properties of to,p or (f> meson which include QCD sum-rule, Nambu-Jona- Lasinio model, Chiral 
perturbation theory or simple hadron effective Lagrangian approach f6V|l0|l . 

In the present work we discuss the vector meson propagation in nuclear matter which could be visualized as meson 
waves traveling through dispersive medium. This concept is usually invoked to study the poinic collective effects in 
nuclear matter Jll[ . A characteristic feature of meson propagation in matter is that there exists two distinct scales 
depending upon its momenta ; if the momenta of the meson is small compared to the Fermi momenta of the nuclear 
matter, the meson propagation becomes inseparable from the macroscopic collective oscillation of the system induced 
by the mesonic field, while on the other hand, for large momenta, as we shall see, the meson propagation can still 
be regarded as usual particle propagation with its dispersion characteristics very similar to that of free space. In 
particular we study the propagation of the uj meson in dense nuclear matter in the long-wave length limit i.e. in the 
regime when the typical characteristic length over which the meson field varies (1/q) is large compared to the the de 



Broglie wave-length of the nucleons 1/pf- Such a study was initially pursued by Chin 1 12 who showed that, in nuclear 
matter, u> meson picks up the collective oscillation of the system and becomes massive. However, in Ref. [ jl2"| , the 
vacuum contribution (Dirac sea) was neglected on the ground that it only renders the couping constant momentum 
dependent and in general, in the long- wave length limit {i.e. for low momentum transfer), the effect is only marginal. 
Later, in a seminal paper by Jean et at JsJ it was shown that, in presence of a mean scalar field (Walecka model), 
the vacuum part contributes substantially with opposite sign. This, therefore, plays a crucial role in determining the 
over all sign and magnitude of the net polarization insertion |jj . In the current work we include both the effects and 
present new and more complete results for the u meson dispersion relations in nuclear matter by considering various 



kinematical limits of the to meson momentum compared to the density characterized by the Debye mass or the plasma 
frequency of the system. Moreover, we also recover the results presented in Ref. [jl2[ in the appropriate limits. 

In addition to the calculation of the dispersion relations, we also examine the analytical properties of the effective 
u! meson propagator for the longitudinal and transverse excitations and construct various sum rules for the relevant 
spectral densities (pt,l)- Later, the full spectral functions (pt.l) are expressed in terms of the residues at the poles 
and cuts of the propagators. Such studies, as we shall see, provide insight to understand and interpret the collective 
excitations of the nuclear matter via the sum rules and the residues in different kinemetic domains of the momenta. 

The paper is organized as follows. We first outline the formalism and describe the collective modes in terms of the 
zeros of the dielectric functions involving density dependent u> meson polarization tensor. In section III, the concept 
of hard nucleon loop (HNL) approximation is introduced which allows for the analytical evaluation of the polarization 
insertions. The dispersion relations for various kinematical domains are then presented. Subsequently, in section IV, 
several sum rules for the lo meson spectral functions are constructed by exploiting the properties of dispersion relation 
techniques of complex variables. Finally the conclusion and summary is presented in section V. 



II. FORMALISM 



It is well known that the collective modes corresponding to density fluctuations are determined from the poles of the 
meson propagators. In Ref. Jl^] it has been demonstrated that in nuclear matter such collective modes are dominated 
solely by the vector interaction in the high-density limit. This, in turn, implies that the in- medium properties of the 
vector mesons are directly linked with the collective excitation set by their propagation. 

The dielectric functions relevant for the to meson propagation in dense nuclear matter can be calculated within 
the frame work of random phase approximation (RPA) which essentially implies repeated insertion of the u> meson 
self-energy involving nucleon- nucleon loop as described in ref. |l2[|T5)| . 

Considering the interaction of the lo meson with the nucleon field, Hi = gytp^f^LO 11 , the second order polarization 



tensor ILj„, can be written as 



n ^ = 7^1 j d "P nH^G(p + q)i^iG(p)} (1) 

where (a, (3) are the isospin indices and G(p) is the in-medium nucleon propagator [T^ ] 

G(p)=G F (p) + G D ( P ) (2) 

where 

Gf(p) = (p^Y + M*)[ *,. ] (3) 
p z — M* z + le 

and 

7 7T 

Gd(p) = dyy" + M*)[— ^(p - E*(p))8( PF - \p\)}- (4) 

The first term in G(p), namely, Gf(p), is the same as the free propagator of a spin 1/2 fermion, while the second 
part, Gd(p), involving 9{pp — \p\), arises from Pauli blocking, describes the modification of the same in the nuclear 
matter at zero temperature. E*(p) = ^/p 2 + M* 2 and M* denotes the in- medium mass of the nucleon, which in the 
present context might be determined from the mean scalar density p3| . 

In a similar vein the polarization insertions can also be written as sum of two parts : 

n^ = n£, + n£, (5) 



From the above set of equations we can write the real part of the density dependent w-meson self energy as |L2 10 

D [ PF d 3 p V^q 2 -Q^{p-q? . . 



^q„)(p v ~ ^q v ) and = {-g^ + ^ 
conforms to the requirement of current conservation, i.e. g^n^ = H? v q v = In order to evaluate n^ 3 conveniently, 



where = (p p — ^#g M )(pjy — ^rqu) and = {—g^ v + ^r 1 )- It is evident that the form of the polarization tensor 



we choose q to be along the z axis i.e. q = (go, 0, 0, | q |), and p.g = E(p)qo— \ p \ \ q \ x, where x is the cosine of the 
angle between k and q. On account of integration over the azimuthal angle, only six components survive. Then again 
for isotropic nuclear matter we have 11^ = 11^ and = 11^. Furthermore, the current conservation condition puts 
additional constraints leaving only two non- vanishing independent component of IT^, , linear combinations of which 
gives us the longitudinal and transverse components of n?„, namely, n£(g) = —Hq + H§ i (q) and (q) = II = 11^ 
The free part is obtained by put ting G (p) = Gf{p) in Eq. ^ which is divergent and can be regularized using 
appropriate renormalization scheme fll2| , |l4f . The regularization procedure we adopt here is the following : 



d n \i F (q 2 )/d{q 2 ) n \ M ^ M ^ =ml = 0, (7) 



where, (n = 0, 1, 2..., oo). This yields 



n^ = 0^« /. to(i-^ og [ M2 _ mMi _ x) ]. (8) 



It is evident that free part does not distinguish between the longitudinal and transverse mode and we can take 
n£ = Ily = H22 = II11. In order to study collective response of the nuclear system and to obtain the dispersion 
relation of the induced oscillation, one needs to find zeros of the dielectric function defined through 

£t,l(9o, \q\) = 1 - — r4 2~nT,L(<7o, \q\) = 0, (9) 

where IIt,l(qo, \q\) — L {q 2 ) + II^ L (qo, \q\)- It is to be noted that the free part Ilf T depends only on the Lorentz 
scalar q 2 while the density dependent part n£ T involves both qo and \q\ individually. The equations CT,L(qo, \q\) = 
in general have to be solved numerically. However, in the next section, we discuss appropriate approximation scheme 
which admits close form solutions. 



III. COLLECTIVE MODES AND DISPERSION RELATIONS 



In this section we discuss first the approximation scheme to derive the polarization tensors analytically. Once the 
polarization insertions are determined one can solve Eq. ^| to find the dispersion relations for the u> meson in nuclear 
matter taking both the effect of free and dense part of the polarization function in the various kinematic limits as 
mentioned in the introduction. In addition, this also provides us a way to interpret the results in terms of the classical 
plasma oscillations. 

First we recall that for collective excitations, the wavelength of the oscillations must be greater than the interparticle 
spacing. This means the meson momenta must be small compared to the nucleon momenta. Therefore, quantitatively, 
it is legitimate to assume that all the loop momenta (nucleon) are hard and the external momenta (meson) are soft 
i.e p ~ pf, \q\ << Pf- We term this as hard nucleon loop (HNL) approximation. This allows us to drop q compared 
to 4 (p. q) 2 in Eq. [| Such an approximation, is, in effect, similar to what has been adopted in Ref. p^ ]. Under this 
assumption the integrations of Eq. [3] relevant for the longitudinal and transverse modes are performed easily to give 
following results : 

to 2 LuLog\^£-} 

and 

Hr(«,*)=2 n + 0- ~ j^t) 2kvp 1 (11) 

where, q = u>, \q\ = k, vp = Pf/^f is the Fermi velocity with ep = \/v 2 p + M 2 and fl 2 is identified as the plasma 

frequency given by £l 2 = fjirgf^. This is related to the Debye screening mass as Q 2 — |m|, [p^ , ^6| . 

Before proceeding further, few comments are in order. It is evident that for k ^ in matter, the longitudinal and 
transverse self-energies are non-degenerate. This gives rise to splitting between these two collective modes in nuclear 
matter which could be attributed to the presence of V^v in Eq. |^. However, they become identical in the static limit 
i.e. n L (w, k -> 0) = n T (w, k -> 0) = Q 2 . 



Next we consider the free part (also known as Dirac sea contribution) given by Eq. Following Ref. ||, this 
can also be simplified further by making a Taylor series expansion of Eq. |8| around M* . To the leading order the 
contribution is found to be P,[l7t 



3tt 2 M 



ni T (w. k) = ^ 77 m 2 + 0(m 2 /4M^) (12) 



As in nuclear matter effective nucleon mass is smaller than its free space value ( M* < M), it is clear from the last 
expression, that the free part reduces the effective ui meson mass in the medium. 

Eq. ^ in conjunction with Eq. O, [l0| and [l^ define the the dispersion characteristics for the collective excitations: 

a , 2 a g* v M*-M 2 3 o2r w 2 , t a; 2 w£ 09 [^±^i 



+ (i - CT ) ;r CTF i = o (is) 



for T mode and 



- - m "-3^^— ^ + 30(l-^ 2 -)[-l+ 2fc ^ ]=0, (14) 

for L mode. 

Generally these transcendental equations, viz. Eq. [l3] and have to be solved numerically, however, analytical 
solutions could be obtained in two limiting cases viz., when k < Q and k > Qby expanding above equations in powers 
of fc and solving the equations by successive iterations. To the first order, when fc << 51 one gets w^r T \ = m* 2 + ft 2 . 
We note that the modes are degenerate in this limit and oscillations are independent of the wave vector k. In the 
case of electron-ion plasma, this mode was identified as the 'plasma waves' or 'Langmuir' oscillation and fl 2 there is 
given by e 2 p/m 2 where p is the density of electron gas and m is the electron mass Jig] . This justifies identifying f2 as 
the plasma frequency in the present pretext. It should be mentioned that here we have absorbed the contribution of 
the free part as 



Once the leading order solution is known, we can solve Eqs. 13 and \\4 iteratively to obtain 



lu 2 = m*J + Q 2 + k 2 + \k 2 v 2 F ."' +•- (16) 

-i = ^ 2 +^ 2 + ^ 2 4^^ + .... d7) 

Evidently, the transverse mode lies above the longitudinal one and it requires more energy for excitation. However, 
the splitting in the low k region is very marginal. In the limit m u — > we recover the results presented in [ p"2"[ , i.e. 

u? T = n 2 + k 2 + -v 2 F k 2 + ... (is) 

5 

uj'I = n 2 + -v 2 F k 2 + (19) 

5 

Furthermore, this taken with Vp —* 1 results in dispersion relations arising out of the photon propagation at finite 
temperature or density (where electron is considered to be ultrarelativistic) with the appropriate definition of the 
plasma frequency |l6| , |l9| ] . It is worthy to note that the contribution of the density dependent part is opposite to that 
of the Dirac part of the polarization tensor. This is consistent with Ref. || . 

Next we consider the case when k > Q but still lu > k and in fact lu ~ kvp. The successive approximation for this 
asymptotic value of k yields solutions of the form 

u T = k 2 + m*J + ^fl 2 + ... (20) 



The corresponding longitudinal frequency takes little bit complicated form; 

2fcV ™* 2 - 2 

3 n 2 1 zn 2 (v 2 F - i 



u£ = k 2 v 2 F (l + 4 e X p{-t^f + ^73^ " 2)) (21) 



Evidently, the longitudinal mode approaches the line lu = kv F and exponentially suppressed for very large values 
of k. We have solved Eq. [l3] and [l4] numerically and checked that for meson momenta k ~ lGeV dispersion results 
are in good agreement with Eqs. II 61 and ITt]. 



IV. SPECTRAL FUNCTIONS AND SUM-RULES 



Equipped with Eq. [To] and [Tl], we can proceed to calculate the spectral functions of the u> meson. As demonstrated 
already, the ui meson in matter has two modes. Confining our attention on the transverse sector first, we write the 
propagator as 

A T (ui,k) = 5 = — — = te ( 22 ) 

w 2_p_ TOi2 _| fi2[ ^ + (1 __^ ) _^ i09[ ^ f]] 

Next chain of steps would be to construct the sum rules by exploiting the analytic properties this effective propagator. 
First we note that the function At is analytic in the complex w-plane having a cut from —kvF to +kvF] in addition 
to the poles ui = ±wt- Having observed this one can make use of Cauchy's theorem for A = At 



A(u,k) 



duj' A(o/, k) 
2ni to' — lu 

x duj' A(oj' + irj, k) — A{u' - irj, k) t / du'A(jJ,k) 

(23) 



2ni u>' — u) /p/ 2iri ui 1 



In the above equation V is a circle with its radius pushed to infinity. This equation can be rewritten in terms of the 
spectral density p(u>, k) containing both the discontinuities across the cuts and the residues at the poles: 

p(uj, k) = 2ImA{u + irj, k) (24) 



as 



A(.,k)=l ^p^fa^A^ 



— oo 



2ir w' — oj .[■£/ 2iri 



It is evident from Eq. that At(u>', k) — > as oj' — > oo, therefore, the integration over the contour T' fails to 
contribute in this region. First sum rule is then immediately obtained by setting lo — in Eq. 

= A T (0,fc) = — — — (26) 



2ir u ' k 2 + to* 2 

Other sum rules are obtained by examining the asymptotic behaviour of At(w,&) when u> — > oo; by taking the 
parity property of the spectral density i.e. pr(w, k) = —pt(—lu, k) into account, one can write from Eq. g5]for uj — > oo 

On the other hand, one can expand Eq. ^2] in powers of w" 1 ; comparing this with Eq. ^ one writes for n = 

UjdbJ 



2?r 



■pr(v,k) = l (28) 



and similarly for n = 1 one observes 



-p T (uj,k) = k 2 +m*J + n 2 (29) 



2tt 



Thus we obtain other two sum rules. Of particular importance is Eq. |28|. This does not depend on k which could 
be shown to be a consequence of the canonical commutation relation. Therefore, it is natural to expect that such a 
relation will also be satisfied by the transverse w meson in nuclear matter. 
The longitudinal propagator is denoted as 



Contrary to the case of transverse propagator, we now have non-vanishing contribution from the contour denoted 
by r' in Eq. as in the limit to — > oo the logarithmic term — > 1. Hence in this limit the longitudinal propagator 
goes like Al(u>, k) — > — 1/fc 2 . Therefore, a subtraction, unlike the transverse propagator, is necessary. The subtracted 
dispersion relation gives the following sum rule essentially in the same manner as in the transverse case. 

cL; p L (u;,k) _ 1 _ 3^ 2 + m*J 

,2tt u " P + M ' j ~ fc 2 (fc 2 + m* 2 + 3fi 2 ) [ ' 

(32) 

Similarly to the asymptotic behaviour of the longitudinal propagator one can write, 

1 i 00 POO J f 1 

This is identical with Eq. [27] except here we have an additional factor of —1/k 2 for reason mentioned in the previous 
paragraph. Likewise, we now can derive the sumrules by first setting n — to get 



/oo 
-oo 



and n — 1 to obtain 



du, , (» 2 +< 2 )(fc 2 +r< 2 + » 2 )-|fc 2 ^ 

p L {u,k) = 2 (35) 



, 2tt 

We can write down the full spectral function of the transverse and longitudinal modes for the uj meson in nuclear 
matter in terms of the poles (and residues at the poles) together with the cuts of the propagator. 

—p T L (uj,k) = Z TtL (k)[S(cu - lo t ,l) - 5{uj + uj t .l)} + Pt,l(v, k)9(k 2 - w 2 ) (36) 

Z7T 



where, /3t,l is given by 



faKfc)= ImIlT - L \A T , L (u;,k)\ 2 (37) 



The imaginary part of the self-energies Ht,l correspond to Landau damping and relevant only for the space-like 
momenta of the meson. Physically this refers to the energy loss due to scattering of the ui meson with hard nucleons 
in the medium. They are given by the imaginary part of Eq.|lO| and EqJTl|. 

The residues at the poles which are used to determine the full spectral function of the ui meson in nuclear matter 
can be determined from the one-loop effective propagators. For instance the residue at the pole for the transverse 
excitation is determined by 

Z T = -QdA^/^U^)])- 1 (38) 

= u T (k 2 v 2 F - c4) 

lo 2 (w 2 - 3(n|, + k 2 + m*J)) + k 2 {k 2 +lj 2 + m*J)v 2 F K ' 

In a similar way we also calculate the reside at the pole for the longitudinal oscillation which is found to be 

z = (fc 2 -^ 2 )(-c|+^fc 2 4) 

L /c 2 [w 2 (-3fc 2 - 3to* 2 +w 2 -3fi 2 ) + fc 2 (fc 2 + m* 2 +w 2 + 3fl 2 )v 2 F ] 1 ' 

We now present the limiting values of these residues at the poles given in section III. First considering k < ft, we 
have for the transverse u 

1 fc 2 5m* 2 + (5 + 34)» 2 

Zr(fc) - 2^ 2 +m: 2 "20 ( (» 2 +TO * 2 )S/ 2 } (41) 

while for large k, i.e. for k > Q we have 



Similar results are obtained for the longitudinal mode also. They are given by, for k < 



iW " 2fc 2 20(tt 2 + m* 2 ) 3 / 2 1 1 

and for k > f2 

ZlW ^ eiK _|^ + g_4_. 2) (44) 

It is now evident from Eq. 28 and Eq. |^| that the transverse spectral function for large meson momentum reduces 
to that of the lo meson propagation in free space, except its mass is modified because of the presence of the mean 
scalar field. Longitudinal mode in this limit decouples from the system. For small momenta, the sum rules are almost 
saturated by the pole contributions. This could be checked numerically. 



V. SUMMARY AND CONCLUSION 



To conclude and summarize, in the present work, we investigate the collective modes induced by the lo meson 
propagation in dense nuclear matter. We then determine the dispersion relations for the lo meson by considering both 
the effect of the Fermi sea and Dirac vacuum within the scheme of HNL approximation. It is observed that in the 
appropriate limit, i.e. for meson momenta small compared to the plasma frequency of the system, the results are 
consistent with the previous calculations |^|,|l2|. The lo meson mass, in presence of a mean field, decreases even at 
normal nuclear matter while for free Fermi gas the effect of nuclcon loop is to increase the lo meson mass proportional 
to the characteristic plasma frequency of the system. 

We also present analytical results for the dispersion relation in the limit k > Q which shows that for the large value 
of the momenta the medium effects die off and the transverse mode propagates like a free meson. The longitudinal 
mode in this limit decouples from the system. We also construct various sum rules satisfied by the lo meson spectral 
density in nuclear matter by investigating the analytic properties of the relevant effective propagator at finite density. 
It is observed that that sum rules are mostly saturated by the pole contributions for small momenta. The residues 
are also presented in the different limit of the meson momenta to understand the collective modes and their behavior. 
Results pertinent to the classical electron plasma oscillations are also recovered. 

Similar studies can also be pursued for other light vector mesons. For instance, it is known that the p nucleon 
interaction is very similar to the lo nucleon interaction, where, in addition to the vector interaction, one has to include 
tensor interaction as well. We, therefore, believe that the p meson propagation in cold nuclear matter can also be 
studied within the framework HNL approximation. This would allow for an analytical derivation for the spectral 
density of the p meson. Hence , sum rules similar to what we present here, relevant for the p meson spectral function 
at finite density can also be constructed in a parallel approach. Such studies are in progress and shall be reported 
elsewhere. 
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